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Abstract 

We present new sets of next-to-leading order fragmentation functions for the production of Kg 
and A particles from the gluon and from each of the quarks, obtained by fitting to all relevant 
data sets from e~^e~ annihilation. The individual light quark flavour fragmentation functions are 
constrained phenomenologically for the first time by including in the data the light quark tagging 
probabilities measured by the OPAL Collaboration. 
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I. INTRODUCTION 



At present, experimental results on inclusive hadron production from e~^e~ collisions 
are the most reliable source for the extraction of universal fragmentation functions (FFs) 
D'^{x,Q'^) (where a labels the fragmenting parton, h labels the produced hadron, x is the 
fraction of the parton's momentum taken by the produced hadron and Q is the factorization 
scale) , which are crucial for making predictions for such processes in future experiments, as 
well as for understanding the non-perturbative mechanism of hadron formation in parton 
jets. However, the extraction of quark flavour separated FFs from experimental data has 
not been completely possible due to the lack of data for processes in which the individual 
light quark flavours are tagged. Therefore, theoretical assumptions between light quark 



flavour FFs had to be made. Recently, we published sets of FFs [1| for each of the three 
light charged hadrons, whose quark flavours were completely phenomenologically separated 
w including the light quark flavour separated measurements from the OPAL collaboration 
2] in the data used for the fltting. This more reliable separation in the light quark flavour 
sector of the FFs via real experimental data is important for the description of hadron 
production in proton-(anti)proton collisions, for example at the RHIC, Tevatron, LHC and 
other experiments, because the proton is composed predominantly of light partons. 

In this paper, we extend our analysis of Ref. Ql to determine FFs for Kg and A pro- 
duction, in which the quark flavours are phenomenologically separated by including, among 
the available data, the tagging probability measurements for each of the quark flavours for 
these two particles provided by OPAL in Ref. P]. FFs for Kg production have been previ- 
ously obtained in Ref. P|, however since no data was available to separate the light quark 
flavours, it was assumed that the d and s quark FFs were equal. Much data for Kg produc- 
tion can also be well described by using FFs for K^ 3, Ql , so it is interesting to verify if 
these two sets of FFs are really consistent. In Ref. P|, FFs for A production were obtained 
by constraining all light quark flavours to be equal and imposing certain relations between 
them and the heavy quark FFs, since no data was available at the time to fully constrain 
the individual quark flavour FFs. A further determination of FFs for A production was 
obtained more recently in Ref. Q], along with FFs for the other octet baryons, again with 
relations imposed between flavours, but also with some assumptions suggested by fermionic 
and bosonic statistics. A more reliable determination of these FFs is important due to the 
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recent data on and A production, taken by the STAR Collaboration j?| at BNL RHIC. 
In particular, next-to-leading order (NLO) calculations deviate considerably from the data 
for A production. For the first time, we present FFs for these particles without imposing 
constraints on the non-perturbative components (except for the choice of parameterization 
for the FFs). 



II. METHOD 



In all cross section calculations in this paper, used to fit FF parameters to data and to 
produce comparisons to this data and other data^we use precisely the same method and 
choice of parameterization, scales etc. as in Ref. and therefore we refer the reader to 
this paper for details. This includes taking the fitted NLO value A^ = 221 MeV of Ref. 

For each parton, our FFs for production are defined to be those for the production 
of a single particle, which is equal to the average of those for K'^ and K^, and our FFs 
for A production are defined to be the sum of those for A° and A°. Our FFs incorporate 
both the intrinsic (non-perturbative) and extrinsic (perturbative, dynamically generated) 
components of hadron production. It is important to note that our FFs contain intrinsic 
transitions involving intermediate hadrons occuring over durations much greater than that 
of the interaction, for example the process g + g^S° + X— i>A + 7 + XinA production 
noted in Ref. j^, since such processes are not subtracted from the data. Treatment of such 
effects are beyond the scope of this work. 

To obtain FFs for production, we fit to all available + e~ — > + X data, covering 

a range of centre-of-mass energies a/s, being the untagged data from TASSO at Jl =14, 
n f — I 

22 and 34 GeVjSi] and at 14.8, 21.5, 34.5, 35 and 42.6 GeV j^], from HRS [10], MARK II 
lul and TPC Q at 29 GeV, from TASSO at 33.3 GeV Q, from CELLO at 35 GeV l3, 
from TOPAZ at 58 GeV Q and from ALEPH Q, DELPHI 1^, OPAL Q and SLD [3 
at 91.2 GeV. The data from DELPHI at 183 and 189 GeV |2o| are not included in the fit 
as will be explained below. In addition, light, c and h quark tagged data from SLD at 91.2 
GeV Il9| and u, ci, s, c and h quark tagged probability measurements from OPAL at 91.2 
GeV \i\ are used to separate the individual quark flavours of the FFs. Likewise, FFs for 
A production are obtained by fltting to all available + A*^(A 
from TASSO at 14, 22 and 34 GeV 181, from HRS U and MARK II 



+ X untagged data 



m at 29 GeV, from 
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TASSO at 33.3 GeV 



ALEPH 



0, DELPHI 2^, OPAL 2^ 



and at 34.8 and 42.1 GeV from CELLO at 35 GeV 



and SLD 



from 



at 9L2 GeV and from DELPHI at 



183 and 189 GeV [20|. The individual quark flavours of the FFs are separated using light, 



9|| and u, d, s, c and b quark tagged 
19|. Soft gluon effects j26| cause the 



c and b quark tagged data from SLD at 91.2 GeV 
probability measurements from OPAL at 91.2 GeV 
DGLAP evolution to fail at small x = 2p/y^, where p is the momentum of the observed 
hadron produced in the final state. Consequently we restrict our analysis to data for which 
X > 0.1, which implies that our FFs will not be valid below this x value. Thus we have 
a total of 193 data points for production and of 129 for A production. The individual 
systematic errors on the data sets are not given in the literature by these experimental 
collaborations, which means that, like all previous FF determinations, we are limited to 
fixing the off-diagonal elements of the covariance matrices for the data sets to zero. We 
stress however that the experimental constraints on FFs (and as{Mz)) would be greatly 
improved by the inclusion of these correlation effects. 



III. RESULTS 

In this section we perform fits to Kg and A production data to determine FFs for these 
particles. In a first analysis, the FF parameters N, a and (3 for the gluon and each of the 
5 quark flavours are released. While an overall good fit is obtained, the fitted initial gluon 
FFs are negative for both particles. Such a result may be unphysical. Redoing the fits with 
the initial gluon FFs fixed to zero (so that the evolved gluon FF in each case is generated 
dynamically from the quarks) gives similar results, implying that the initial gluon FFs are 
consistent with zero with respect to these data. This is to be expected since the gluon FF 
only enters the cross section at NLO, and can therefore only be reasonably constrained by 
more accurate data such as that in Ref . ||li] . To obtain a set of FFs in which the gluon FF is 
reliable, we redo both fits with the initial gluon FF for each particle fixed to some function 
whose choice and motivation will now be explained. Using the approximation 

= iff. (1) 

where b = u,d if a = d,u, otherwise b = a, we fix the inital gluon FF for Kg production to 
half the AKK jlj| initial gluon FF for K~^ + K~ production. For A production, it is reasonable 
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to assume that the gluon FF is related to that of the proton. However, since the valence 
structure of A is uds compared to uud for the proton, we assume that this relation holds 
with an overall suppression factor, which we determine approximately by requiring a good 
description of the pp initiated cross section measurements from the STAR collaboration [3] 
to be discussed later. We find that a good description of such data is achieved if the initial 
gluon FF for A production is fixed in the fit to the AKK initial gluon FF for p+p production 
suppressed by a factor of 3. Thus the STAR data is in some sense crudely included in the 
fit, at least in the sense of providing some constraint on the gluon FF for A production. 

We obtain = 1.14 and x^f = 1-39 for the fits to the Kg and A production data, 
respectively, indicating that both fits are good overall. The quality of the fit to each individ- 
ual data set is determined by its Xdf value, and these are listed in Tables HI -II VI The values 
in Tables H] and |n] show that the agreement with the OPAL tagging probabilities is good. 
The good agreement with heavy quark tagged data is surprising when compared with our 
findings of Ref . . For that analysis of light charged hadron production data, although we 
found good agreement with the DELPHI, SLD and TPC heavy quark tagged data, a poor 
description of the OPAL heavy quark tagged data was obtained, which, as we suggested, 
may be due to large angle gluon emission effects. Whatever effect caused this disagreement 
is clearly not as significant for the K'g and A production data we are considering here, al- 
though we note that we find a rather large Xdf value for h quark tagged production. 
The Xdf values for production in Table IIIII indicate that each data set is well fitted 
except the DELPHI data at 183 and 189 GeV which are excluded from the fit. When these 
data are included, their x^p values remain high (and not much less than those shown in 
Table HIT]) . Such data require further consideration. In any case, only 3 data points of these 
DELPHI data lie in the region x > 0.1, so they do not significantly change the size of the 
data sample used in our fit. We will therefore exclude them in our fit. Table HVl shows that 
unsatisfactory agreement was obtained only with untagged data from DELPHI and b quark 
tagged data from SLD, both at = 91.2 GeV. The data points from DELPHI at 183 GeV 
(1 point) and 189 GeV (1 point) for A production were excluded for the same reasons as in 
the case of Kg production, although the point at 183 GeV could be well fitted. 

The results for the FF parameters are shown in Table IVl The value for a and/or f3 may 
become large because the shape of the FF at small and/or large x respectively is not well 
constrained. However, the value of is then large to compensate for this. For A production, 
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the inital FF for the d quark is negative, probably because this initial FF is consistent with 
zero with respect to the data which constrains it. The initial u quark FF has a large negative 
a, due to a lack of data to constrain this FF at small x. 

The graphical comparison with the quark tagging probabilites is shown in Fig. The 
curves obtained from our FF set in this analysis (labelled "AKK" ) agree well with the data. 
For Kg production, we also show the curves obtained from the FF set for production 
in Ref. assuming the relation in Eq. (^. These curves are consistent with these mea- 
surements. Also shown are the curves obtained from the FF set for Ka production of Ref. 

n 

|27| (labelled "GR"). These curves are consistent with the data, showing the validity of 
these authors' method of obtaining these FFs from the Monte Carlo generator HERWIG. 
We note that curves calculated from the FFs for production from Ref. (also ob- 
tained from HERWIG) using the same assumptions (Eq. (0)) are also consistent with the 
Kg quark tagging probabilites. For A production, we also show the curves obtained from 
the FF set presented in Ref. (labelled "FSV"). These are generally consistent with the 
data except for the s and c quark tagged data. Unfortunately, we are unable to show the 
predictions from the FF set for A production of Ref. Q, since the Mellin transform of the 
parameterization used therein, that is required for our Mellin inversion calculation of the 
cross section {]J|, cannot be obtained in closed form. Such predictions require a NLO x space 
program. The comparison with the differential cross sections is shown in Figs. |2] and El For 
Kg production, the disagreement of the calculation with the DELPHI data point at 183 
GeV can be seen, and, of the two points at 189 GeV, with the one that is the higher in x. 
For A production, the failure of the description of the data point at 189 GeV can be seen. 
Note that the DELPHI data for A production at 91.2 GeV that gave a high Xbf rather 
small errors. Finally, in Fig.Elthe cross section for A production in the region x < 0.1 is not 
shown since large negative values were obtained there. This results from the large negative 
value for a in the FF for the u quark discussed earlier. 

The quark tagged data from SLD is shown in Figs. ID and El together with the calculation 
of these cross sections from the FF sets of this analysis. Reasonable agreement is found with 
most of the data. However, in Fig. El there is significant disagreement with the interval from 
the b quark tagged cross section for A production whose central point is just below x = 0.4. 
The light flavour quark tagged cross section in this figure is large and negative in the region 
X < 0.1, again due to the large negative value for a in the FF for the u quark. 
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We find that the large initial strange and charm quark FFs for A production render the 
other quark FFs negligible at large x, and this is reflected at higher scale in the OPAL tagging 
probabilities at large Xp. This finding was not observed in Ref. ^, where the light quark FFs 
were constrained to be equal. This should make a significant difference to hadron production 
data from proton- (anti)proton collisions, where the light quark PDFs are important. 

To test our FFs, we use them to predict the recent preliminary measurements for pp ini- 
tiated cross sections at ^/s = 200 GeV from the STAR Collaboration 7^ and for pp initiated 
cross sections at ^/s = 630 GeV for Kg and A production from the UAl Collaboration 12811 
at the CERN SPS. We use the NLO coefficient functions for this process from Ref. (29l |. 
and the CTEQ6M parton distribution functions from Ref. [s^. Figure IHl shows the results 
for Kg production, which are very similar to those obtained in a similar plot in Ref. as 
is to be expected from the similarity between the two AKK curves of Fig. ^ While good 
agreement with STAR data is found, some other effects are required for the description of 
the UAl data. Figure [3 shows the results for A production. As discussed earlier, the STAR 
data were used to motivate the choice of the initial gluon FF: Firstly, fixing the initial gluon 
FF to that of the AKK proton from Ref. [l| and performing a fit to the e'^e~ data, we found 
that the description of the STAR data was too high by a factor of about 3. Since we also 
found that these data are dominated by the gluon FF at low factorization scale, we divided 
this choice of the initial gluon FF by this factor and redid the fit (i.e. fitted the initial quark 
FFs) to get our final FF set for A production. The plot shows good agreement with both 
STAR and UAl data within the theoretical errors. 



IV. CONCLUSIONS 

We have obtained FFs for Kg and A production by fitting to data from e+e^ collisions 
over a wide range of ^/s values, from 14 to 91.2 GeV. To separate the light quark flavour 



FFs, and to improve the determination of t 



le heavy ones, we have included for the flrst time 



the quark tagging probabilities from Ref. 



with the FFs from Ref. 



• 3- 



For Kg production, we flnd good agreement 



and the FFs for K^ + K^ production from Ref. |l| (using Eq. 
in the latter case) adding support to the reliability of the data of Ref. P]. For A production, 
we found that the light quark flavour FFs are signiflcantly different at large x, with the 
initial strange quark FF dominating over the remaining quark FFs. 
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In both cases the initial gluon FF is fixed in the fit. For production, we use the AKK 
initial gluon FF for (divided by 2) . This leads to a good description of the pp initiated 
STAR data for i^Spro duct ion, but a bad description of the UAl data. This problem was 
also found in Ref. l| using the AKK FFs for production via Eq. For A production, 
we use the AKK initial gluon FF for p/p production, after dividing by a factor of 3. This 
factor is chosen to give a good description of the STAR data for A production, and also leads 
to a good description of the UAl data. Since such data constrains the gluon well at low 
Pt, it should in future analyses be directly included in the list of data fitted to. However, 
currently the calculation, which implements Monte Carlo integration, is too slow to obtain 
fitted FFs in a reasonable time. 

In order to make predictions, our fitted FF's over the range 0.1 < x < 1 
and Mq < Mf < 200 GeV can be obtained from the FORTRAN routines at 
http://www.desy.de/~simon/AKK2005FF.html, which are calculated using cubic spline in- 
terpolation from a linear grid in x and linear interpolation from a linear grid in In y/s. 
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TABLE I: Xtif values obtained from the measured quark tagging probabilities for K'^ production, 
ryf , at ^ = 91.2 GeV in Ref. ||. 



a 


d 


u 


s 


c 


b 


Xdf 


0.27 


1.38 


0.54 


1.77 


2.36 



TABLE IL Xdf values obtained from the measured quark tagging probabilities for A production, 
r?^, at ^ = 91.2 GeV in Ref. \ 



a 


d 


u 


s 


c 


b 


Xdf 


1.05 


1.00 


0.29 


1.39 


1.02 



10 



TABLE III: CM energies, types of data and Xdf values for various samples of Kq production 
mesurements. Samples not used in the fits are marked by asterisks. The columns are labelled by 
the quarks that were tagged in the mesurements {q implies no tagging). 



tagged: 

[GeV] 


q 


uds 


c 


6 


14.0 


1.82 [2] 








14.8 


2.65 [9j 








21.5 


0.69 [9j 








22.0 


0.99 [2] 








29.0 


1.27 0.40 0.16 [12] 








33.0 


0.82 [13] 








34.0 


1.33 [8] 

mm 








34.5 


1.89 [9] 








35.0 


0.30 f9] 1.64 [14] 








42.6 


1.02 [2] 








58.0 


0.01 [1^ 








91.2 


0.33 0.57 [17] 0.31 [18] 1.13 [19] 


0.93 [19] 


0.54 [19] 


1.85 [19] 


183.0 


21.87 [22]* 








189.0 


4.41 [20]* 
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TABLE IV: CM energies, types of data and Xdf values for various samples of A production 
mesurements. Samples not used in the fits are marked by asterisks. The columns are labelled 
by the quarks that were tagged in the mesurements {q implies no tagging). 



tagged: 

Vi [GeV] 




uds 


c 


b 


14.0 


0.73 [8j 








22.0 


0.73 [2] 








29.0 


1.19 [21] 0.83 [22] 








33.0 


1.59 [13] 








34.0 


1.85 [2] 








34.8 


2.29 [23] 








35.0 


0.88 [14] 








42.1 


0.90 [22] 








91.2 


0.43 4.05 [24] 0.47 [25] 0.33 [19] 


1.63 [19] 


1.54 [19] 


3.70 [19] 


183.0 


4.89 [20]* 








189.0 


16.34 [20]* 









12 



TABLE V: Values and errors of A^, a and /3 resulting from the fit. 



Hadron 


Flavour 


N 


a 


(3 




d 


0.0297 


-0.949 


0953 




Vij 


1.29 


0.137 


4.53 




s 


0.560 


-0.283 


1.43 




c 


8.43 


0.550 


6.04 




b 


28.8 


0.574 


13.5 




g (fixed) 


7.96 


2.72 


2.45 


A 


d 


-16.7 


-0.169 


20.2 




u 


0.0016 


-3.67 


4.27 




s 


46043 


9.27 


9.11 




c 


156 


2.71 


9.51 




b 


244 


2.15 


15.6 




9 (fixed) 


0.289 


0.130 


0.854 
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AKK 

AKK[(K"+K")/2, u<->d] 

GR 

FSV 



FIG. 1: Quark tagging probabilities rj^ixp, s) at ^/s = 91.2 GeV. The solid curves labelled "AKK" 
are calculated from our FF set for production in the graphs on the left hand side and our 
FF set for A production in the graphs on the right hand side. The curves of the same label from 
the FF set in Ref. ji] (after interchanging u and d) are also shown. In addition, we show the 
prediction for Kg production from the Kg FF set of Ref. 27] (labelled "GR" ) and the prediction 
for A production from the A FF set of Ref. jsf (labelled "FSV"). The corresponding measured 
OPAL probabilites of Ref. || are also shown. 
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FIG. 2: Normalized differential cross section of inclusive Kg production. The curves are calculated 
from the FFs obtained in our analysis, at the various energies of the data shown. Each curve and 
the corresponding data are rescaled relative to the nearest upper one by a factor of 1/5. 
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FIG. 3: Normalized differential cross section of inclusive A production. The curves are calculated 

from the FFs obtained in our analysis, at the various energies of the data shown. Each curve and 
the corresponding data are rescaled relative to the nearest upper one by a factor of 1/5. 
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FIG. 4: Normalized quark tagged differential cross sections for inclusive Kg production compared 
with corresponding data from SLD The curves are calculated from the FFs obtained in our 
analysis. Each curve and the corresponding data are rescaled relative to the nearest upper one by 
a factor of 1/20. 




FIG. 5: Normalized tagged differential cross sections for inclusive A production compared with 

□ 

corresponding data from SLD [19||. The curves are calculated from the FFs obtained in our analysis. 
Each curve and the corresponding data are rescaled relative to the nearest upper one by a factor 
of 1/20. 
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FIG. 6: The invariant differential cross section for inclusive Kg production in pp and pp collisions. 
Preliminary data from the STAR Collaboration 7] at ^/s = 200 GeV (divided by a factor of 30 
for clarity) and data from the UAl Collaboration 28] at ^/s = 630 GeV are shown, together with 
their predictions using the FFs obtained in this paper. In the latter case, the upper, central and 
lower curves are calculated with a renormalization and factorization scale of ^ = Mf = pt/"^, Pt 
and 2pT respectively. 
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